Dispersive stabilization of the inverse cascade for the Kolmogorov flow 
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It is shown by perturbation techniques and numerical simulations that the inverse cascade of kink- 
antikink annihilations, characteristic of the Kolmogorov flow in the slightly supercritical Reynolds 
number regime, is halted by the dispersive action of Rossby waves in the /3-plane approximation. For 
/3 — * 0, the largest excited scale is oc lnl/f3 and differs strongly from what is predicted by standard 
dimensional phenomenology which ignores depletion of nonlinearity. 
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Planetary-scale flow is subject to the competing effects 
of quasi-two-dimensional turbulence and Rossby waves 
(/3-effect). It is known from phenomenological arguments 
and numerical simulations that the inverse cascade which 
characterizes the large-scale dynamics of two-dimensional 
turbulence in planetary flow can be halted by Rossby 
wave dispersion and that the ensuing flow exhibits alter- 
nating jets. jl|-|5||- The standard argument of Rhines 
rests on a comparison between the local eddy turnover 
time and the period of Rossby waves. More generally, 
the interaction of waves and turbulence is a subject with 
a wealth of applications in astro/geophysical flow and 
plasmas (see, e.g., Ref. ||). 

The case of strongly dispersive waves is reasonably well 
understood through the theory of resonant wave interac- 
tions (Rcfs. |6[0] and references therein). It is however 
inapplicable to a rotating planet where the dispersive ac- 
tion of Rossby waves (with a frequency (3/k) is felt only 
at the very largest scales (small wavenumber k). For the 
case of weakly dispersive waves, the nonlinear dynamics 
in the absence of waves must be understood in detail be- 
fore we can find how they are affected by the presence 
of the waves. We cannot resort just to dimensional ar- 
guments which ignore the depletion of nonlinearities, an 
effect which is very common in turbulence |§ and which 
can take extreme forms in some situations, as we shall 
see. Our present study deals with the Kolmogorov flow 
H for which the inverse cascade is understood rather 
well in terms of kink-antikink dynamics [To[| . The basic 
Kolmogorov flow is u = (cosy, 0), obtained by applying 
to the two-dimensional Navier-Stokes equation a force 
f = ^(cosy, 0). The Kolmogorov flow has a large-scale 
negative eddy viscosity instability when the kinematic 
viscosity v is just below the critical value v c = l/v2 
P). The large-scale dynamics are then, to leading order, 
one-dimensional QJllJl- Close to v c and in the presence 
of Rossby waves, the large-scale dynamics are governed 
by the one-dimensional "/3-Cahn-Hilliard" equation with 
cubic non linearity, derived by multiscale techniques in 
Ref. |7j (see also Ref |ll| for the case (3 = 0), which reads 



Here, U(v) = s 2 (v 4 / (2T 2 ) — v 2 ) is a quartic potential and 
d^ 1 denotes spatial integration for zero- average functions 
in the interval [0,L] over which periodicity is assumed. 
(In the original setup of Ref. M, the constants are s = 
1/V3, T = i/3/2 and A 3 = 3/V2.) 
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FIG. 1. Simulation of the Cahn-Hilliard equation with 
13 = 0. (a) : evolution of the total energy, (b) : evolution 
of the the energies of various Fourier modes (as labeled). An 
inverse cascade is observed from iV = 8 to TV = 3. Eventually, 
the dominant mode becomes TV = 1 (not shown in the figure) . 

We begin with qualitative considerations illustrated by 
simulations. For /3 = 0, the solutions to this equation live 
essentially within a slow manifold of soliton-like solutions 
with an alternation of plateaus at the zeros v = ±r of the 
potential, separated by alternating kinks and antikinks 
|l2| . The corresponding fixed points, having TV pairs of 
regularly spaced kinks and antikinks, are all unstable sad- 
dle points of a Lyapunov functional, except for TV = 1 
which gives a stable absolute minimum. The temporal 
evolution is a cascade of annihilations of kink-antikink 
pairs, leading eventually to the gravest TV = 1 mode [ fl3| . 
This is illustrated in Figs. 0and||(a) obtained by numer- 
ical simulation using the method described in Ref. 0. 
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In the Fourier space an inverse cascade is observed with 
the dominantly excited wavenumber shifting roughly to 
smaller and smaller wavenumbers (see, e.g. Ref |Q and 
Fig. |l|(b)). Except for the short kink-antikink annihila- 
tion episodes, the motion of kinks is described by simple 
ODE's to be derived below. They involve exponential 
couplings between adjacent kinks, so that the typical du- 
ration of the plateaus increases exponentially as N de- 
creases. 
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FIG. 2. Solutions with four kink-antikink pairs, (a) : snap- 
shot of velocity v(x, i) in a slowly evolving state for f3 = 0; 
(b) : final stable asymptotic state for f3 = 10 -3 with TV = 4 
exhibiting distorted plateaus between the kinks (this is actu- 
ally a traveling wave measured in a suitable moving frame). 
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FIG. 3. Same as Fig. [j] but for (3 = 10~ 3 . Propagating 
Rossby waves accelerate the cascade and, eventually, halt it 
when it reaches N = 4 (four pairs of kink-antikinks) . 

For f3 0, the dispersive action of the waves mod- 
ifies the cascade. For small values of /3, the solution 



retains the characters of kink dynamics with superim- 
posed propagating waves which are expected to predom- 
inantly modify the dynamics at large scales where the 
Rossby wave period is the smallest. As the destabilizing 
coupling between the kinks decreases exponentially with 
the typical separation of adjacent kinks, the stabilizing 
effect of dispersive waves can dominate for solutions of 
sufficiently large wavelength, thereby halting the inverse 
cascade at an intermediate wavelength A = L/N with 
integer N > 1. This is illustrated with a simulation for 
f3 = 1CP 3 shown in Fig. |3| where the cascade is found 
to halt at wavenumber k = 4, corresponding to a wave- 
length A = L/4. 

We turn now to more systematic theory, beginning 
again with (3 = 0, case for which we mostly follow 
Kawasaki and Ohta JrJ . The (pure) Cahn-Hilliard equa- 
tion admits M{x) = ±rtanhsx as stationary solutions 
of d%Mj — U'(Mj) = 0. When the typical separation 
between kinks is large, a solution with N kink-antikink 
pairs may be written, near the jth kink, as v(x,t) = 
Mj(x) + Vj(x,t), where Mj{x) = ejM(x — Xj) with 
€j = (— X)- 1 and exponentially small remainder iij(x,t). 
Hence, the time derivative is, up to exponentially small 
terms, dtv(x,t) = — Y^e=o 1 ie(t)d x Me(x). Substituting 
into (Q), integrating twice in space, multiplying by d x Mj 
and integrating over the periodicity interval [0, L], we ob- 
tain 

— y~] xt j / d x Mj(x)g 2 {x - x')d x M e (x')dxdx' = 

^3 £=0 Jo JO 

^ ^v](3Mj+v j )d x M j dx-2e j rh(t), (2) 

where h(t) is a time-dependent integration constant 
which can be determined by the constraint of momen- 
tum conservation. Carrying out the integrations in (g), 
we obtain the ODE's for kink motion: 

Ajixt = e-'to-"'- 1 ) - e-sOw-ay). (3) 

Here, Xj is the location of the j-th kink and Aje is a 
symmetric matrix given by 

^ = 8A^ { ej -« " Xl) + lk?) - h Sj - e ) ' 

where Q2(x) is the L-periodic Green's function satisfying 
d x Q-z(x) = —S(x). When L — > 00, Qijxj — xt) reduces to 

— \xj — xe\/2 and A is easily inverted p3|. Although the 
dynamics described by (||) borrows most of its character 
from the heteroclinic connection between unstable fixed 
points, the trajectory in phase-space does not generally 
proceed from the vicinity of one fixed point to another. 
Fast jumps from one slow manifold to the next may occur 
at distance from the fixed points. Thus, as the number 
of kinks and antikinks decreases with time, the dominat- 
ing wavenumber does not necessarily decrease monoton- 
ically as happens in the "arithmetic" inverse cascade of 
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Rcf. [[14j. This may be seen in Fig. |l|(b) where the dom- 
inant mode goes from wavenumber k = 8, to 5 and then 
to 6. 

The stability of a given fixed point w (0) of @j 
with respect to kink motion is obtained by differen- 
tiating (H). It is convenient to define the Fourier 
components ip m of kink displacements 5xj as ip m = 
(2N)- 1 J2 2 J^S 1 Sx 3 e^ (m ^ N \ After some algebra, we 
obtain ip m = croVVn for m < N, with the eigenvalue <jq 
given by 
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with 9 m = nm/N and A = L/N (the subscript refers 
to (3 = 0). The corresponding eigenvectors for a given 
m are v a (x) = Y^j=o C0S j9md x M(x — Xj) and 

v b{x) = Y^q 1 i- 1 ) 3 sin j8 m d x M(x - Xj). The leading 
order of (|j) corrects a factor 2 error in Ref. jl0| in which 
the exponential variation of Vj near Xj+i and Xj—i was 
neglected. Note that our result is consistent with Ref. 
and agrees with our simulations. For large A, all the 
eigenvalues <to (actually cto(to)) for m > are positive, 
thereby demonstrating the instability of multi-kink-pair 
solutions. 

We now turn to the case of non-vanishing small (3, 
which is studied by a singular perturbation method that 
we outline. The stationary solution goes over into a 
(slowly) traveling wave solution which, in a suitable 
frame moving with the velocity c = (3c\ + (3 2 c 2 + 0((3 3 ^ 



may be written in the time-independent form v 
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(3v (1) + (3 2 v {2) + 0{(3 3 ). Here, T> (1) and v (2 *> satisfy 
FtfV = Q {a) and Fv^ = Q« with T = d 2 - U ', 

g<°) = A 3 - x (cxd- 1 ^ - a- 3 #°)) , and Q« = 

A3- 1 (dfl-v 1 ) - a-w> + c 2 o-^) + \k (« (1) ) 2 , 

where Uq = U(v^). At first order, c\ is obtained from 
the solvability condition J d x v^ dx = as 
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where A — 2.404113- • -. Like ctq in (jij), c\ is given in 
(0) up to exponentially small terms (for large A = L/N) 
and has excellent agreement with numerical simulations. 
(Note that, at next order, C2 vanishes.) The first-order 
perturbation v^- 1 ' of the traveling wave profile can also 
be obtained perturbatively in a large-A expansion (not 
given here). This leads to a distorsion of the plateaus, as 
seen in Fig. ||(b). 

The equation governing the perturbation Sv(x, t) = 
v(x + ct, t) — v(x) can be written as 



d t (p = £ip, ip= d x 1 Sv, 

£ = ~X 3 d x (d 2 - U"(v))d x + cd x - /3d- 1 . 



Unlike the case with (3 = 0, the perturbation to the 
stationary solution of the /3-CH equation does not re- 
duce simply to a change in the kink locations. We 
must also consider the dispersive effect of the /3-term 
which modifies the shape of the slow modes and con- 
tributes to stability. Substitution in (||) of the expan- 
sions of c and v in powers of [3 leads to an expansion 
£ = Co + (3£i + (3 2 £ 2 + 0{(3 3 ). The unperturbed eigen- 
value (To becomes a = ctq + i(3fii + (3 2 a 2 + i(3 2 [i 2 + 0((3 3 ). 
Similarly, the eigenfunctions <p a — d x ~ 1 v a and ipb = d~ 1 Vb 
are respectively modified as tp a + (3<p a i + [3 2 (p a2 + 0{(3 3 ) 
and ipt, + (3<pt,i + (3 2 (pb20([3 3 ). A hierarchy of equations 
is then obtained for ip ai ipb, fail Vfcii l Pa2i ■■■■ As is 
usual in such singular perturbation problems, the cor- 
rections to the eigenvalues are obtained from solvability 
conditions. This gives (with (/, g) = L^ 1 J Q L fgdx) 

((fbXlCPa) = -Hl{lph,iPb), 

(tpa,£l<Pal) + (<Pa,£2<Pa) = Cr 2 ((p a ,(p a ) - fXl{<p a , fbl) • (?) 

The coefficient /ii gives only the shift in the imaginary 
part of the eigenvalue. For the real part, the second-order 
term a 2 is needed. It is obtained from (Q); limiting our- 
selves to leading order large-A contributions, we obtain 
after considerable algebra 
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Note that the correction is negative and, hence, stabiliz- 
ing. Though the effect is small, it increases algebraically 
with A, while the nonlinear coupling of kinks decreases 
exponentially in (Q). Therefore, stabilization of the m- 
mode perturbation to the stationary solution is obtained 
at leading order for large A when a + a 2 [3 2 < 0, that is 



> 0c = J 35,389,440- 
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The condition is the most restrictive for m = 1, that is 
q = tan7r/(2A). 

When A is only moderately large, an accurate esti- 
mate of (3 C requires, in principle, a careful determina- 
tion of (pbi and ifiai by matching of solutions near to and 
away from the kinks. We have used an alternative semi- 
numerical approach, which works for all values of A, in 
which the auxiliary equations stemming from the pertur- 
bation theory, are solved numerically with a discretiza- 
tion over a quarter- wavelength of ip a and ipb, using a 
procedure written with Mathematica, available from the 
first author (BL) . The eigenvalues are obtained from the 
solvability conditions of the discretized problem, where 
degeneracies are lifted by imposing symmetries with ap- 
propriate boundary conditions for each linear problem to 
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be solved. We have also solved directly the stability prob- 
lem by differentiating a Galerkin-truncated expansion of 
(|l|) in Fourier modes. Newton's algorithm is used to ob- 
tain the traveling solution and its phase velocity; eigen- 
modes being then calculated with a QR algorithm. The 
latter procedure provides both the fast and slow modes 
with exponential separation of the eigenvalues. As the 
problem becomes very stiff when L is large, we are in 
practice limited to values of L less than a few hundreds. 

Table I compares the asymptotic and the semi- 
numerical solutions of the perturbation theory with the 
direct solution of the stability problem. There is excel- 
lent agreement between the direct calculation and the 
semi-numerical solution of the perturbation problem for 
N < 5. It breaks down when the approximation of sep- 
arated kinks is no longer valid. The asymptotic result 
( |To| ) provides then the correct order of magnitude but is 
wrong by a factor 2 or 3. The convergence to the nu- 
merical solution is observed at larger L but is very slow : 
multiplying L by 10 and 50 narrows the discrepancy for 
mode 2 to 4.6% and 1.3%, respectively. 

Numerical temporal integrations of ([!]) starting from 
random initial conditions eventually achieves a traveling 
wave solution with one of the values of N found stable 
by the perturbation theory. 

The main result of our perturbation theory is ( |To| ) 
which gives the minimum value of able to stabilize 
a solution of period A = L/N with N > 1. By in- 
verting A in terms of /3, we can explain the halting of 
the inverse cascade at a wavelength A which scales as 
A = — (2/s) In (to leading order). Note that stan- 
dard phenomenology, based on dimensional analysis a 
la Rhines Q with equilibration of nonlinear and Rossby 
characteristic times, gives a drastically different scaling, 
namely A ~ |/3| -3 . (The nonlinear time is oc A 2 and the 
Rossby time to l/(/3A); velocities are 0(1).) The dis- 
crepancy arises from the failure of dimensional analysis 
to capture the almost complete suppression of nonlincar- 
ity obtained in the plateaus. 
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-35.002 2.358 10~ b 


2.358 10" b 


-35.002 1.45 10" b 
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-16.067 4.965 10" 4 


4.966 10~ 4 


-16.067 2.20 10" 4 
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-9.2073 1.069 10 -2 


1.067 10~ 2 


-9.2090 3.2 10 -3 


5 


-5.9507 8.96 10 -2 


10.28 10~ 2 


-5.9637 1.8 10~ 2 



TABLE I. Comparison of the perturbation theory in both 
its asymptotic and numerical form with the direct stability 
calculation. The calculation is done for L = 76.93 (10 unsta- 
ble modes when f3 = 0) and several values of N . 



The one-dimensional character of the large-scale dy- 
namics of the Kolmogorov flow in the slightly super- 
critical regime amplifies the discrepancy. In multi- 
dimensional high-Reynolds numbers turbulence, coher- 
ent structures (vortices, filaments, sheets, . . . ) exhibit 
also strongly depleted nonlinearities which dimensional 
arguments fail to capture. It would be of interest to 
study dispersive stabilization for a strongly nonlinear 
high-Reynolds number inverse cascade of the kind con- 
sidered by Kraichnan [Q. So far no systematic theory 
can handle this, although some progress has been made 
recently on a related passive scalar problem [jL7j . 

Finally, we mention that in recent work, devoted to 
zonal jets in planetary flow, Manfroi and Young |l§| l 
studied a closely related problem in which stabilization 
is provided by a friction term — rv added to the right 
hand side of ([!]). We find that a suitable adaptation of 
the analysis for the dispersive case gives stabilization for 
r > r c ~ e~ sA /A. 

We are grateful to Joanne Deval for careful checking 
of all the calculations and to an anonymous referee for 
useful remarks. 
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